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Thin film YIG (Y 3 Fe 50 i 2 ) is a prototypical material for experiments on thermally generated pure 
spin currents and the spin Seebeck effect. The 3ct; method is an established technique to measure 
the cross-plane thermal conductance of thin films, but can not be used in YIG/GGG (Ga 3 Gd 50 i 2 ) 
systems in its standard form. We use two-dimensional modeling of heat transport and introduce a 
technique based on Bayesian statistics to evaluate measurement data taken from the 3ct; method. 

Our analysis method allows us to study materials systems that have not been accessible with the 
conventionally used 3ct; analysis. Temperature dependent thermal conductance data of thin film 
YIG are of major importance for experiments in the field of spin-caloritronics. Here we show data 
between room temperature and 10 K for films covering a wide thickness range as well as the magnetic 
field effect on the thermal conductance between 10 K and 50 K. 


Knowledge of the thermal properties of a system is fun¬ 
damental in experiments including heat flow. In solids, 
heat is mainly transported by long-wavelength acoustic 
phonons. Therefore, even comparably thick ’thin film’ 
systems encounter the generic issue that the heat con¬ 
ductivity depends on the film thickness, since long wave¬ 
length phonons will be scattered by the surfaces or in¬ 
terfaces. For the electrical conductance such an effect 
is described by the Fuchs-Sondheimer relation.EI^ The 
derivation of this relation was simplified by the compa¬ 
rably easy access to experimental data, i.e. measurement 
of electrically conducting thin films on insulating sub¬ 
strates. In addition, due to the Pauli principle only elec¬ 
trons at the Fermi surface participate in transport, so 
that a Fermi velocity and mean free path for the the¬ 
oretical analysis can be computed. For heat transport 
the full phonon spectrum needs to be considered and the 
heat conduction of the substrate can usually not be elim¬ 
inated in the experiment rendering the determination of 
thin film heat conductance an ambitious task. 

The interaction of the spin degree of freedom with 
an energy flow in a solid, spin-caloritronics, is a re¬ 
cent field of researcfP^ in which YIG (Y 3 Fe 5 0i2) is 
of interest due to its low magnetic damping. Temper¬ 
ature driven pure spin currents and the spin Seebeck 
effect have been observed using YIG based sample^^l^^ 
and there is an ongoing discussion about the origin 
of the spin Seebeck effect,!^ which includes magnon- 
phonon interaction. Here the so far unknown tempera¬ 
ture dependence of the thermal conductance of thin film 
YIG is needed to quantitatively determine the real tem¬ 
perature gradients in the investigated systems. While 
values of the bulk thermal conductivity vary between 
(6 ... 8) W/(m Padture et al.t^ discuss that their 

result of approx. 6 W/(mK) possibly underestimates the 
true value by up to 1.5 W/(mK). Theoretical calcula¬ 
tions of the magnon thermal conductivity in bulk YIG 
exist ,11^ but vastly disagree in magnitude with experi¬ 
mental result s.l^ 


Mostly for numerical reasons conventional methods to 
determine thin film thermal conductances fail in sys¬ 
tems in which the film thermal conductivity is compa¬ 
rable to that of the substrate. Therefore, crucial temper¬ 
ature gradients a re us ually calculated using bulk ther¬ 
mal conductivitie^^^ even if films of several hundred 
nm thickness are used. In this paper we introduce a 
data evaluation scheme based on Bayesian statistics to 
determine the thickness dependent heat conductance of 
thin film samples. Our novel method allows us to de¬ 
termine the thermal conductivity of systems not only in 
which the substrate has a heat conductivity much larger 
than that of the film material, but independent of this 
condition. We present the thickness dependent thermal 
conductivity of thin film YIG and demonstrate that at 
temperatures far below room temperature the thermal 
conductance is mainly limited by point defects. 

To determine the thermal conductance we applied the 
‘5uj method that uses a patterned metallic structure as 
a heater and meas ureme nt device for the thermal resis¬ 
tance of the sample .11^^ An alternating current with fre¬ 
quency uj is applied to the heater structure, heating the 
sample and modulating the resistance R with a frequency 
of 2uj. Ohm’s law and trigonometric identities imply a ‘5uj 
contribution to the voltage, [/s^, related to the temper¬ 
ature oscillation in the heater. furthermore depends 
on the heater’s temperature coefficient of resistance and 
the electric heating power applied. Generally, the ther¬ 
mal penetration (probing) depth of the experiment de¬ 
pends on the thermal diffusivity of the material and the 
applied frequency as , so that with increasing fre¬ 
quency the thickness of the layer ‘seen’ decreases. While 
at low frequencies the experiment is sensitive to both 
film and substrate, at high frequencies the thermal pen¬ 
etration depth decreases until only the film is measured. 
The thermal conductivity can then be inferred from the 
slope of with respect to In (2cc;) (slope method) .1^^ 
Fig. [T] shows the temperature oscillation in dependence 
of In (2uj). In the high-frequency range above 50 kHz the 
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FIG. 1. The temperature oscillation in the heater structure in 
dependence of In (2ct;) measured between 100 Hz to 100 kHz 
with cj the dimensionless angular frequency measured in Hz. 
The inset shows the temperature dependence of the resistance 
of the An heater strip (black) and its temperature derivative 
(red). 


data of the 6.7 /im and 2.1 jam films overlap, indicating 
that the thermal penetration depth is smaller than the 
thickness of the thinner film (at 50,525 Hz the thermal 
penetration depth amounts to 2.1 /im). Particularly in 
this regime the finite heater width needs to be taken into 
accounfl^. In the standard procedure at low frequencies 
the heat resistance of the film is compared to that of a ref¬ 
erence sample with a smaller film thickness (differential 
method). This approach requires the thermal transport 
to be strictly perpendicular to the sample surface by as¬ 
suming that the thermal conductance of the substrate by 
far exceeds that of the film and that the thermal pene¬ 
tration depth is larger than the film thickness. The first 
requirement concerns the choice of the substrate. The 
second can only be met by an adequate heater structure 
and measurement frequency. 

For YIG the choice of substrates is largely limited to 
other garnets. However, the traditional 3uj method is not 
usable for YIG films on Gd 3 Ga 50 i 2 (GGG) substrates, 
as the thermal conductivities of film and substrate are 
comparable .1^ To conserve the validity of the classical 3uj 
approach, the information depth of the method needs to 
be smaller than the film thickness, which requires elec¬ 
tron beam lithography and high-frequency detection (see 
Fig. [^. In addition, for very high frequencies and very 
low temperatures the amplitude of the temperature os¬ 
cillations gets very small and the slope based method 
becomes inapplicable (below approx. 50 K in the case 
of YIG). Thus, at sufficiently low temperatures not only 
do the assumptions of the differential 3uj method, e.g. 
the substrate as a heat sink, break down, but also the 
applicability of the slope method. This underlines the 
necessity for a new approach to determine thermal con¬ 
ductances. 

The film’s thermal conductance and the temperature 


oscillation in the heater are connected by the solution of 
the heat diffusion equation. Instead of remaining in the 
strict cross-plane (one-dimensional) transport regime, 
the requirements for the heater geometry and the sub¬ 
strate thermal conductivity can be relaxed by extending 
the equation for the temperature increase in the heater 
to a model of two-dimensional (in- and cross-plane) heat 
transport.l^^ll^For an arbitrary system of films on a sub¬ 
strate, the solution for the complex temperature oscilla¬ 
tion in a heater line is given by the integral equatiorP^ 
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where P is the applied heating power and the heater 
stripe has the dimensions 2h and 1. Ai and Bi are given 
by recursive sequences Ai and Bi depending on the in¬ 
plane and cross-plane thermal conductivities tvxy,i and 
f<iz^i^ the thickness di of the z-th layer of the stack and 
the respective thermal diffusivities ai = given by 

PiCs,i 

the thermal conductance, the density p and the specific 
heat Cg. As arbitrary thicknesses di can enter the cal¬ 
culation, this method is independent of the thickness of 
the individual layers. However, as the thermal ‘wave’ 
decays exponentially within a layer, layers thicker than 
approx. 100 pm are difficult to analyze. Interface ther¬ 
mal (Kapitza) resistanc es are n eglected here, but can in 
principle be introducedj ^^l"^ l "^ ^ 

The real part of the temperature oscillation resem¬ 
bles an improved version of the traditional approach 
by including the anisotropy of the thermal conductance. 
Therefore, in principle it is possible to determine also 
the in-plane contribution to thermal transport without 
a second 3uj structure and without varying the width of 
the heater line. Furthermore, the real part depends on 
the density and the specific heat only by cross-terms in¬ 
troduced by the multiplication of the complex quantities 
Ai and Bi. The imaginary part in particular includes an 
a~\ term and therefore is a measure for the density and 
the specific heat, but is also sensitive to the cross-plane 
thermal conductance. 

If for each layer {tixy^ p^ Cg} are unknown, the equa¬ 
tion contains 4 • Yiayers free parameters and it could the¬ 
oretically possess a unique solution, if the measured set 
of frequencies was equal to or greater than the number of 
parameters. The large number of parameters then leads 
to a high degree of computational complexity and insta¬ 
bility. Assuming all properties of all layers except for 
that of interest are known from previous measurements, 
the number of parameters can be reduced to those of the 
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film. A sensitivity analysis of Eq. to the material pa¬ 
rameters allows for further simplification. We found that 
any determination of p and will be prone to large un¬ 
certainties and is only possible for materials with very 
large density and specific heat capacity. Thus, determin¬ 
ing the thermal conductivities of the film can be simpli¬ 
fied by assuming reasonable values for the film’s density 
and temperature dependent specific heat capacit^ttSl and 
only determining hixy and t<iz. 

Fitting Eq. [^to measurement data is complex due to 
the integral equation. Therefore, techniques s uch as a 
differential material properties search algor it hnPSEH or 
neural network^^ are implemented to find a meaning¬ 
ful fit. It may also be possible to use the Monte Carlo 
technique to obtain a solution, which, however, does not 
need to be unique. E.g. the algorithm presented by Ol¬ 
son et al.l^ lacks stability and the results depend on the 
initial conditions of the fit. To circumvent this difficulty, 
we have implemented a new algorithm based on Bayes’ 
theorem, Eq. Bayesian statistics uses previously deter¬ 
mined information (‘prior probability distribution’) and 
then uses data to compute a ‘posterior distribution’. 



FIG. 2. Temperature dependent thermal conductances of 
YIG samples (first batch: 6.7 pva (red circles), 2.1 /xm (black 
squares) and 190 nm (blue triangles)) and exponential fits to 
the high-temperature regions (solid lines). The inset displays 
the fit by the Gallaway model around the maximum. 
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Here rui denotes the Tth model of a series of models, 
o a measurement and p{mi\o) the conditional probabil¬ 
ity that a model rrii describes the physical reality pro¬ 
vided the observation o. Vice versa, p(o|m^) is the non- 
normalized likelihood that an observation o be made pro¬ 
vided rrii is real. The denominator is the effective nor¬ 
malization of the probability p{mi\o) to [ 0 , 1 ]. Assuming 
an experiment is carried out and the outcome is com¬ 
pared to Njn models of the process, as a prior one could 
choose a uniform distribution. The difference between 
model value rrii and observational value o can be ex¬ 
pressed in terms of the uncertainty of the observation a 
by = I I. This number represents an approxima¬ 
tion for the probability that the observed value actually 
does not come from that model, so that 


p{o\mi) = 2 - 20(Ei), (6) 

where 0(E) is the cumulative Gaussian distribution run¬ 
ning from 0.5 to 1 , as the E^ are positively defined. 

For each combination of the parameters {rixy^ 
within a physically meaningful range the expected AT 
is calculated to form an extensive catalog of models for 
temperatures between room temperature and 10 K. As 
a naive prior a uniform probability distribution function 
(PDF) for hixy and riz is chosen and then improved by dis¬ 
favoring a large anisotropy in the thermal conduct ance.lSI 
Each experimental data point at a specific frequency is 
compared to all models at that frequency using Eq. 
and Eq. This leads to an ‘intermediate’ PDF, which 
is recursively used as a prior for the next frequency data 


set and thereby improves the total PDF with increasing 
amount of data. The PDF determined after all data has 
been included represents the posterior PDF for n^y and 
i<iz from which the expectation value and the variance are 
computed. 

To determine the credibility of our analysis method, we 
re-evaluated data taken on Lao.eTCao.ssMnOsP^, which 
had been evaluated using the conventional differential Scj 
method. Our results agree very well within the experi¬ 
mental uncertainty and reproduce the results obtained 
in the cited reference. In addition, we compared a sam¬ 
ple consisting of Au/Al 2 03 /Mg 0 to a second sample 
Au/MgO to determine possible interface contributions. 
The results reproduce literature values of the thermal 
conductivity of AI 2 O 3 without the need for an additional 
thermal interface resistance and therefore Kapitza resis¬ 
tances can be neglected. 

We employed our analysis method to YIG grown on 
(5x10) mm GGG substrates. A first batch of sam¬ 
ples was deposited by liquid phase epitaxy (LPE) and 
pulsed laser deposition (PLD). The film thicknesses were 
6.7 pvn and 2.1 /xm (LPE) and 190 nm (PLD). The second 
batch included samples grown by LPE with thicknesses 
of 1.51 /xm, 3.08 /xm, 7.78 /xm, 12.01 /xm, 22.83 /xm and 
50.66 /xm. The gold layer from which the heater struc¬ 
ture was fabricated was grown by DC sputter deposition. 
AFM measurements show the YIG film roughness to be 
below 1 nm. X-ray diffraction spectra demonstrate that 
sputtered An has a polycrystalline fee structure and X- 
ray reflectivity results in a RMS roughness of 4 nm. The 
heater strip (thickness: 50 nm) is rectangular in shape 
with a length of 1 mm between the voltage contacts and 
a width of 20 /xm, which yields a room temperature spe¬ 
cific resistivity of An of 2.6 • 10“^ Om. The temperature 
dependence of its resistance and temperature coefficient 
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of resistance are shown in the inset of Fig. and fol¬ 
low the Bloch-Griineisen law. In particular, the slope of 
the resistance with the temperature is a measure of of 
the signal resolution, so that below 10 K no meaning¬ 
ful 3uj measurement can be performed. Caused by the 
low temperature coefficient below 10 K a measurement 
of the thermal conductivity is only possible above this 
temperature. An Anfatec eLockIn 205/2 lock-in ampli¬ 
fier, assisted by a Wheatstone bridge and common-mode 
subtractioiP^, was used to apply the heating voltage and 
read out 1 / 3 ^. 

A slope-based measurement of the thermal conduc¬ 
tance of a GGG substrate was used to supply the al¬ 
gorithm with the thermal conductance of the substrate. 
While the slope method yielded a cross-plane thermal 
conductance for GGG of i<iz = 7.8 W/(mK), the two- 
dimensional method resulted in nz = (8.6±1.6) W/(mK). 
The room temperature thermal conductance of YIG was 
determined for two samples from the high-frequency 
data, in which the reduced thermal penetration depth 
causes the method to be sensitive to the film only. The 
cross-plane thermal conductance of the films is thereby 
determined to f<iz = (9.0T0.2) W/(mK) at room tempera¬ 
ture. The Bayesian fitting algorithm was then applied to 
the low-frequency data below 2 kHz and led to values of 
i^xy = (9.5T1.1) W/(mK) and Kz = (8.5T0.6) W/(mK), 
which demonstrates the reliability of our algorithm. 

The Bayesian approach is capable of providing the tem¬ 
perature dependent thermal conductance down to low 
temperatures. Fig. shows the thermal conductance of 
the first batch of YIG films and displays clear peaks. 
At high temperatures, Umklapp scattering is the dom¬ 
inant factor that reduces the thermal conductance and 
the thermal conductances are expected to be thickness 
independent, while the peak is a measure for the length 
scale at which phonon boundary scattering becomes pre¬ 
dominant. Therefore, a smaller film thickness should lead 
to a shift of the peak towards higher temperatures and, 
correspondingly, to lower peak values of the thermal con¬ 
ductance, which is observed in the first batch of samples. 
In bulk material Boona et al.^ found a peak thermal con¬ 
ductivity of 110 W/(mK) at 25 K. At high temperatures 
the data fit an exponential decay close to the maximum 
as is expected from phonon transport theory. 

Modeling of the low-temperature behavior requires an 
in depth analysis of contributing transport processes. 
While we could experimentally determine the temper¬ 
ature dependence of the thermal conductivity, a detailed 
analysis of the different scattering contributions will re¬ 
quire an independent and extensive determination of the 
defect structures for each sample, which is outside the 
scope of the current paper. Nevertheless, we have per¬ 
formed a phenomenological deconvolution of different 
scattering contributions in order to point out the chal¬ 
lenges and possibilities. Specifically, thermal conductiv¬ 
ity data can be modeled by contributions from different 
relaxation times r as described by Callawa}^^. In par¬ 
ticular, scattering of phonons at frequency uj at point 
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FIG. 3. Peak temperature Tmax (upper plot) of the thermal 
couductauce aud peak value A^max (lower plot) as fuuctiou of 
film thickuess. The data poiut at 500 /xm is takeu from Booua 
et al.^ The black squares result from the first batch of sam¬ 
ples, while the red circles correspoud the secoud batch of sam¬ 
ples. The error bars are derived from the sample-to-sample 
variatiou of Tmax and /^max- The blue aud red hues iu the 
lower plot iudicate the depeudeuce ou the defect couceutra- 
tiou parameter of the maximum iu the thermal couductivity 
at coustaut scatteriug parameters other thau the thickuess 
effect. 


defects is described by = Auj^ with A a constant 
proportional to the defect density. The samples of the 
first batch grown by LPE can be explained using a single 
set of relaxation times and only varying the thickness- 
dependent boundary scattering relaxation time based on 
the actual thickness of the films (lower panel of Fig. [^. 
Data taken from the film grown by PLD deviates from 
this model. A different parameter set is required in which 
the largest difference between the samples is a different 
defect concentration given by the model parameter A. 
The order of magnitude difference in the values of A im¬ 
plies a large difference between the defect densities. The 
defects do not need to be of the same type in both sets of 
samples and a difference in film quality (PLD resulting 
in a larger defect concentration) is to be expected. In the 






5 



FIG. 4. Magnetic field dependent contribution of the thermal 
conductivity of YIG films with thicknesses of 50.66 fim and 
1.51 fim. 

case of the bulk sample data taken from literature there 
may be additional influences that are unknown to us. 

It can be expected that films approaching the bulk 
state approach bulk values of k, but this effect is unob¬ 
served in our data (Fig.[^. There is no apparent system¬ 
atic trend, even though all thermal conductivity curves 
can be fitted individually using the Callaway model, but 
with strongly varying defect concentrations. The sample- 
to-sample variation determined from two randomly cho¬ 
sen samples within one wafer of the maximum thermal 
conductivity amounts to 8% and that of the maximum 
temperature to 6%. This confirms the reproducibility of 
our results, as the 3uj method generally is prone to sys¬ 
tematic uncertainties of up to 10%. Our data imply that 
the quality of the film, i.e. the defect concentration, is a 
vital piece of information required to understand the thin 
film thermal conductance of YIG. At room temperature 
the values of the thermal conductance obtained using our 
model are slightly larger than the literature value. The 
fact that the slope-based values of the thermal conductiv¬ 
ity of YIG and the result of the presented algorithm agree 
well at room temperature implies that our new method 
is capable of providing physically sound results. 

Finally, Fig. [^displays the difference in thermal con¬ 
ductivity him of YIG films between 10 K and 50 K in 
samples with thicknesses of 50.66 fim and 1.51 /im when 


in a magnetic field of 8 T, which freezes out the magnonic 
contribution,f^and in zero field. Statistical errors are de¬ 
termined by sampling 100 data points per temperature 
point, which is necessary due to the small magnitude of 
the difference. Our data show that Km may be thick¬ 
ness dependent and that the magnetic field may affect 
the thermal conductivity differently in thick films, while 
the thin film does not show a contribution significantly 
different from zero. Above 40 K the difference cannot 
be distinguished from thermal noise. The data suggest 
that the temperature of the largest difference is at ap¬ 
prox. 15 K in the thick film, while there is no significant 
difference in the thin film at all. The point of largest 
difference is at a temperature much lower than the peak 
thermal conductivity, which is in agreement with Boona 
and Heremans^^. Generally, we infer that the magnetic 
field dependent contribution to the thermal conductiv¬ 
ity below 50 K does not exceed 1 W/(mK). Further re¬ 
search will be required to eliminate systematic uncertain¬ 
ties, which led to small but unphysical negative values of 
the thermal conductivity, and to understand the mag¬ 
netic field effect on the thermal conductivity in YIG. 

In conclusion we developed a method determining the 
thermal conductance of thin film samples based on two- 
dimensional heat transport and Bayesian statistics. We 
applied the method to YIG films on GGG substrates 
where the original 3uj method is inapplicable. For thin 
film YIG we determined a room temperature thermal 
conductance of Kz = (8.5 ± 0.6) W/(mK) with low 
anisotropy. We present temperature dependent thermal 
conductances of thin film YIG between room tempera¬ 
ture and 10 K with sample thicknesses between 190 nm 
and 50 fim. The results agree with bulk values at room 
temperature, but differ significantly towards lower tem¬ 
perature, where the thermal conductance is shown to be 
significantly reduced at small thicknesses compared to 
bulk material. In particular, in films grown by LPE the 
effect of the defect density exceeds that of the expected 
thickness dependence. 
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